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RECIPROCITY RELATIONS IN AERODYNAMICS 1 


By Max. A. Heaslet and John R. Spbeiteb 


SUMMARY 

Reverse-flow theorems in aerodynamics are shown to he based 
on the same general concepts involved in many reciprocity the- 
orems in the physical sciences. Reciprocal theorems for both 
steady and unsteady motion are found as a logical consequence 
of this approach. No restrictions on wing plan form or flight 
Mach number are made beyond those required in linearized 
compressible-flow analysis. A number of examples are listed, 
including general integral theorems for lifting, rolling, and 
pitching wings and for wings in nonuniform dovmwash fields. 
Correspondence is also established between the buildup of circu- 
lation with time of a wing starting impulsively from rest and 
the buildup of lift of the same wing moving in the reverse 
direction into a sharp-edged gust. 

INTRODUCTION 

Some of the most important results V. the recent study 
of wing theory have been achieved through the development 
of reverse-flow relations. The theorems already obtained 
are of outstanding practical utility and it appears obvious 
that the fullest exploitation of the methods has yet to be 
accomplished, either from a purely theoretical standpoint 
or in the routine calculation of wing characteristics. Atten- 
tion to such problems in aerodynamics was initiated by von 
Kdrman (ref. 1) who first announced the invariance of drag 
with forward and reversed directions of flight for a nonlifting 
symmetrical wing at supersonic speed. Subsequently, 
advances in the theory were made by Munk, Hayes, Brown, 
Harmon, and Flax (refs. 2 through 7). Up to the present 
time, the most general results have been expressed by 
Ursell and Ward (ref. 8) and by R. T. Jones (ref. 9) in his 
attack on the study of wing shapes of minimum drag. 
In the forms given in the two latter papers, the derived 
equivalences could be termed reciprocal or reciprocity rela- 
tions rather than reverse-flow relations ; in fact, this change 
in terminology divorces attention, momentarily, from the 
purely aerodynamic aspects of the results and, in this way, 
suggests a reorientation in terms of the various similar 
relations appearing in other engineering fields. In the 
theory of elasticity, for example, a reciprocity theorem for 
small displacements of an elastic medium is so expressed 
as to appear in formal agreement with the statement of the 
result given by Ursell and Ward (see, e. g., ref. 10). This 
theorem is attributed to E. Betti and was published in 1872. 
A generalization was given by Lord Rayleigh in 1873, and 


in various sections of his two volumes on the theory of sound 
(ref. 11) discussions of reciprocal relations in an elastic 
medium and for acoustic sources are given. In 1886, von 
Helmholtz (ref. 12) obtained, by means of variational 
methods applied to Hamilton’s characteristic function, a 
reciprocal theorem for small changes in the momenta and 
coordinates of a general dynamical system in forward and 
reversed motion. This result was commented on, in turn, 
by Lamb (ref. 13) and an independent proof based upon 
Lagrange’s generalized equations was given. The paper by 
Lamb is of particular interest since it contains the essential 
idea underlying the development of reverse-flow theorems 
in wing theory. Thus, Lamb remarks, as had Lord Rayleigh 
previously, that reciprocity relations between sound sources 
do not apply directly in a moving atmosphere. He points 
out, however, that the reciprocity can be restored if the 
direction of the wind is also reversed. Further examples 
of reciprocal theorems appear in the theories of electricity 
and magnetism (in particular, reference should be made 
to Maxwell’s discussion of the subject in ref. 14) and of optics. 

The generality in the statement of reciprocity relations 
appears, almost universally, to have held back their applica- 
tion to problems for which they are obviously, in retrospect, 
particularly fitted. This generality is even more apparent 
in some of the conclusions of Lord Rayleigh and von Helm- 
holtz which apply to nonconservative systems. 

The purpose of the present paper is twofold. First, a close 
connection will be established between reverse-flow theorems 
in subsonic and supersonic, steady-state wing theory and 
known reciprocity relations between two solutions of the 
equation governing the flow .field. In this way, machinery 
will be provided whereby extensions of existing results to the 
case of unsteady motion follow directly. Second, a number 
of particular problems in wing theory in steady and unsteady 
flow will be considered. It will be shown that, provided 
attention is limited to force and moment characteristics, 
the complexity of many solutions involving nonuniform • 
flow fields, control-surface deflections, and unsteady motion 
can be reduced considerably. In some cases, previously 
obtained solutions will be calculated. Comparison with 
the original calculations will almost invariably highlight 
the economy of effort in obtaining the final result. The 
utility of reverse-flow theorems is based on the fact that 
they build from known solutions and thus avoid the necessity 
of .starting each problem anew. 


1 Supersedes NAOA TN 2700, “Reciprocity Relations In Aerodynamics” by Max. A. Heaslet and John R. Sprelter, 1052. 
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GENERAL ANALYSIS 


RECIPROCITY RELATIONS FOR A CLASS OF PARTIAL 
DIFFERENTIAL EQUATIONS 


, In this section, integral relations associated with linear 
partial differential equations will be reviewed from the stand- 
point of relating independent solutions. The subject 
matter is precisely Green’s theorem and, in common with 
the usual expression of the theorem, it is preferable to 
treat the variables initially as abstract quantities. Consider, 
therefore, a class of linear partial differential equations of 
second order with independent variables X u X 3 , . . . , X„ 
that may be thought of as rectangular coordinates in a 
space of m dimensions. Denote differentiation of the 
function yp(X u X 3 , ... , X m ) with respect to the variables 
X t and X] by the subscript notation 

M^bxjbXj 

and consider the differential equation 


L(if) aSS A t] (yp) t] +Byp= 0 (1) 

f=i J=i 


where, for the purposes at hand, A {j —A Jf are made inde- 
pendent of X t and X h and B is a constant. Such equations 
fall within the class of self-adjoint equations. 

By Green’s theorem (see, e. g., ref. 15), it is possible to 
relate two arbitrary functions yp and G by the integral expres- 
sion 


JJ/WW) - ypL(Q)]dV= _JJ [QD n p- ypD n Q] dS (2) 


where the left member is a volume integral over a prescribed 
region in m-dimensional space and the right member is a 
surface integral over the hypersurface S enclosing the given 
region. Equation (2) certainly holds for any region in 
which yp and Q and their first and second derivatives are 
continuous. The directional derivaties D n are defined in 
terms of the direction cosines n it . . ., n m of the normal to 
the surface S with the stipulation that the normal is directed 
into the given region. Thus, setting 


X) njAi^Nvi 


where v u v s , . . ., v m are the direction cosines of a line termed 
the “conormol,” the directional derivative is defined by the 
axpression 

(3) 

i-1 ° v 

If, finally, \p and O are assumed to satisfy equation (1), the 
left side of equation (2) vanishes and the resulting expression 


to the particular case when yp and 0 are identified ■with the 
perturbation velocity potential <p in the theory of incom- 
pressible-fluid flow. The governing equation is Laplace’s 
equation in three dimensions 

= 0 ( 6 ) 

and the reciprocity relation takes the form, for two possible 
solutions <p and <p', 

where the directional derivatives are now along the true 
normals to the surface S enclosing a three-dimensional 
volume. It is known that any actual state of motion of a 
liquid for which a single-valued velocity potential exists can 
be produced instantaneously by the application of a properly 
chosen system of impulsive pressures. These impulsive 
pressures are directly proportional to the velocity potential 
plus an arbitrary constant which may, in the present case, 
be associated with the pressure of the uniform stream. 
Equation (6) is thus seen to represent summations of the 
cross products of impulse and normal velocity in two possible 
motions of a conservative system and is a special case of the 
dynamical theorem (ref. 11, p. 98) 

(7) 

1 r-l 

where p r , q T and ])/ , g/ are generalized components of 
impulse and velocity in any two possible motions of a system 
which starts from rest. 

The interpretation of equation (6) that leads to equation 
(7) provides an indication of the close connection between 
reciprocal theorems based upon the principles of least action 
and the symmetric character of Green’s theorem for certain 
second-order differential equations. In the subsequent 
applications it will be convenient to proceed directly from 
equation (4) and seek to establish reciprocal relations between 
flow fields in wing theory. Such a process is well known 
when yp or G is replaced by the elementary solution associated 
with a unit source and, in this case, establishes a general 
solution in terms of source and doublet distributions deter- 
mined by arbitrary boundary conditions. The present 
objective is, however, different in that one wishes to get a 
symmetrical dependence between two general solutions. 
Moreover, the apparent symmetry of equation (4) must be 
consistent with physical considerations so that, for example, 
the velocity at point A induced by a singlo source at point B 
is equal to the velocity at point B induced by a source at 
point A. If, in the two systems, the effective free streams or 
flight directions are opposed, a fore and aft symmetry occurs 
and the possibility of maintaining symmetry in the reciprocal 
relations becomes feasible. 


^ ^QJD n ypdS=^ ^ypD n Q.dS 


(4) 


is a general reciprocity relation expressing the functional 
dependence between two arbitrary solutions of equation (1). 

An interesting interpretation of equation (4) has been 
given by several writers (see, e. g., ref. 16, p. 46) and applies 


RECIPROCAL THEOREMS IN WING THEORY 

REVERSE FLOW FOR SUBSONIC WINGS 

Consider a thin wing at possibly a small angle of attack and 
situated in the immediate vicinity of a flat surface Miich is 
designated the xy plane. For sufficiently small thickness and 
angle of attack, the perturbation velocity potential or any of 
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tho perturbation velocity components of the wing satisfy the 
linearized partial differential equation of compressible flow. 

/SVa+l/W+^K— 0 (8) 

where /3 2 = 1 —M 0 3 = 1 — ( Uo/a^) 1 , and U 0 ,a-o are, respectively, 
the flow velocity and speed of sound in the free stream. 
Equation (8) applies in forward or reverse flow, provided the 
corresponding free-stream Mach numbers M 0 and M 0 ' axe 
equal. In figure (1) a lifting wing with plan form P is 
indicated along with the vortex wakes as the configuration 
would appear if the two flow fields were superimposed. It 
will be assumed that the wing chord is finite and that the 
profiles are closed. 


and 

ip dS=^ \j/^dS; n directed downward (9b) 

<ri <ti 

where the primes denote conditions for the reversed flow. 

In equations (9), \p is now replaced by the perturbation 
velocity potential <p(x,y,z) and by the x or streamwise 
component of perturbation velocity u(x,y, z). By virtue of 
the irrotationality of the flow, the gradients of u and w are 
related by the expression 

bu_ bw 
bn bx 




It is proposed to apply the general reciprocity relation (4) 
to these flow fields in a manner similar to the approach used 
in developing the basic solutions of the differential equation 
(see, e. g., refs. 17 and 18). Thus, for subsonic flow, hemi- 
spherical regions of large radius and lying first above and 
then below the plane of the wing will be chosen as the 
volumes of integration. The surface integrals will therefore 
extend over a hemisphere with center at the wing and a flat 
surface that lies immediately adjacent to the z = 0 plane. 
This latter surface is subsequently to be brought into coin- 
cidence with the plane of the wing but must be considered 
first in its displaced position since only then can the flow be 
assumed free of possible singularities in perturbation veloci- 
ties and their gradients. As in wing theory, in general, the 
attenuation of the perturbation potential and its gradients 
may be assumed of such a nature that the integrated con- 
tributions of tho wing and its wake over the hemispherical 
surfaces vanish in the limit as the radius becomes infinitely 
largo. It remains, therefore, to consider the integrals over 
the surfaces at z= 0+e and z—O—e. Denote these surfaces 
as <r h and <n, respectively, where the subscripts u and 2 specify 
values above and below the 2=0 plane. Equation (4) can 
then be written as 



n directed upward 


(9a) 


and equations (9) can both be written in the form 
Integration by parts over either of the surfaces <r„ or a t gives 

f* 00 

At x =— co the potential <p for the forward flow vanishes and 
at z= + c» the upwash w' in the reversed flow vanishes so 
that the first term on the left is zero. , The remaining double 
integrals have for their surfaces of integration the displaced 
planes <r u and <rj. In order to obtain a concise form of the 
reverse-flow theorem, it suffices to subtract the integrals ex- 
tended over, crj from the integrals over <r u and let the planes 
approach coincidence with the plane of the wing. Since 
w,vf and u,u' are continuous everywhere except possibly 
in the immediate vicinity of the wing, the integration region 
can be restricted to planes slightly above and below the wing 
but extending beyond the wing edges. Provided the singu- 
larities at the edges can be disregarded, the analytical ex- 
pression of the steady-state reverse-flow theorem of Ursell 
and Ward (ref. 8) becomes 

(Uu'Wu' — UfiDi)dS=j j (Uu'w u — Ui'lOi)dS . (11) 

P / 

for either lifting surfaces or symmetrical wings where P is the 
plan form of the wing in the 2=0 plane. 

It remains to discuss the effect of edge singularities. In 
the case of a lifting surface, square-root singularities in both 
u and w can occur at the leading and trailing edges just on 
and off the wing, respectively. In the combined flow fields, 
the limi ting process would then yield residue terms analogous 
to the leading-edge thrust of a lifting -plate. If, however, 
the Kutta condition is imposed at the trailing edge for both 
the flow in forward and reverse direction, a combination of 
singularities does not occur and equation (11) is valid. If 
the leading edge of a symmetrical wing is rounded so as to 
produce a square-root singularity inwon the wing, a square- 
root singularity in u occurs just off the wing and a term cor- 
responding to the leading-edge drag (ref. 19) appears. If 
the geometry of the wing is fixed in the forward and reverse 
flow, however, the effect of these terms is canceled and equa- 



• dy — J* fi uw'dS=ij ^u'wdS (10) 
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tion (11) still applies. It is important to point out that if 
the Kutta condition is not satisfied, then the area of integra- 
tion of equation (11) cannot be restricted to the plan form. 

The two sides of equation (11) are expressed in terms of the 
same coordinate system but it is usually preferable to asso- 
ciate with each of the two streams an x axis extending in the 
stream direction. To this end introduce now the subscripts 
1 and 2 to denote the forward and reverse flow and the two 
coordinate systems. Thus, in general, 

5Ci=— Xa+f, Vi=— Vt+V, Zl=3l ( 12 ) 

where £ and ij are arbitrary constants, and equation (11) then 
becomes 


the subscripts 1 and 2 refer to forward and reverse flow in 
the two axial systems. 

, REVERSE FLOW FOR SUPERSONIC WINGS 

The development of a reverse-flow theorem for supersonic 
wings parallels closely the analysis for the subsonic case. 
For either planar or multiplanar systems, the conormal in 
equation (4) is, in fact, the normal so long as the surface of 
integration is a plane parallel to the x axis. In the case of 
the single wing, for example, equations (9) apply where the 
surfaces a are slightly removed from the plane of the wing 
and where ip satisfies the differential equation 

(Mo 1 — 1) ypa— 'Pn— = 0 


Pi 


J J [«u t (3> v y^) (*i,2/i) ~'>ii 1 (xi,y\)wi a (xuyi)}dx l dy l = 

JJ ['Muj fayi) W U1 (Xi,y s ) —U h (x 2 ,y») w,j (z^j)] dx^dy* (13) 
In the case of a symmetrical nonliftin^ wing, the relations 


u u =u h w u ——w t 

must apply and, in the case of a lifting surface, linearized 
theory yields 

U u ——Ui, Wu—Wt 

It follows that in either case, equation (13) reduces to the 
form 


JJ «i (xi,yi)w 3 (zi,i/i)<iridyi== u? (x 2 ,y 2 ) w i fe,^) dxdv* (14) 
Pi Pi ■ 


where the velocity components can be evaluated on either 
the upper or lower surface of the plan form. If, further- 
more, the linearized pressure relation 

p — Pa— — poUqU ’ (15) 

is used where p is local static pressure, p 0 is static pressure 
in the free stream, and the wing profiles are assumed closed, 
equation (14) becomes 1 


JJ (x u yi)wi(xi,yOdxidyi=JJ p 2 (a^y,) (x?,y 2 ) dxdy 2 (16) 

Pi Pi 


If, instead of specifying boundary Conditions in a single 
plane, it is necessary to treat boundary conditions for a 
system of planes, the expression of the reverse-flow theorem 
is of the same general form as equation (11). Provided the 
Kutta condition is imposed at the trailing edges of all lifting 
surfaces, the relation becomes 

JJ (F„WS=JJ (V n ) 2 UidS (17) 

2 2 


where the area of integration 2 extends over both sides of 
all the wing surfaces, V n is the component of perturbation 
velocity normal to and directed away from each wing, and 


In the limit as a approaches the s=0 plane, the reversibility 
theorem takes the form of equations (11) and (16), provided 
the integration extends beyond the edges of the wing. It is 
necessary to include these edges for wings with subsonic 
leading and trailing edges since singularities occur in the 
perturbation components and . the solutions are not neces- ' 
sarily unique. For supersonic-type edges, the area of 
integration can be confined to the plan form of the wing and 
this is also true for subsonic edges, provided the Kutta con- 
dition holds for all subsonic trailing edges in both the forward 
and reverse flow. Equation (17) relates the two possible 
flows in the case of multiplanar systems. 

REVERSE UNSTEADY MOTION 

In the case of unsteady motion at either subsonic or super- 
sonic flight speed, the basic equation may be taken in the 
form 

— = 0 (18) 

where t—Oxf/, do is the speed of sound in the undisturbed 
region of the field, f is time, and 4> is the perturbation velocity 
potential or any of the perturbation velocity components. 
Equation (18) is the acoustic equation for small disturbances 
in three space dimensions and' holds for a system of Cartesian 
coordinates fixed relative to the undisturbed air. In appli- 
cations to wing theory, therefore, the wings move relative to 
fixed axes. 

In the derivation of a useful theorem it is convenient to 
treat thin wings at small angles of attack and to assume that 
the motion takes place in the xy plane. The visualization of 
the time and geometry relations is relatively easy for two- 
dimensional wings moving at a uniform speed, as indicated in 
figure 2. The airfoil starts at time t'=0 and moves to the left 
at a constant velocity U 0 so that the trace of the leading edge 
in the xt plane is x= — Uot'= —Mo f, and the trailing-edge trace 
is x=2a—Mot. The lines z= ±t and x=2 a±t are the traces 
of the extremeties of the regions affected by the acoustic 
waves set in motion at t'= 0 by the leading and trailing edges. 
In figure 2, the wing has traveled a time T'=T/a 0 and the 
boundary condition determining the ■wing shape during the 
motion will be fixed by prescribing the value of vertical 
induced velocity w over the region “swept out” of the xt 
plane by the w ing - In order to determine a reciprocal theo- 
rem, a second wing is assumed to start at the final position of 
the first wing and to move with negative velocity until it has 
reached the initial position of the first wing. With these 
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Fioore 2. — Two-dimensional wing In it plane. 

concepts in mind, it follows from equation (4) that the 
relation 

\l5t wdxdt ( 1Q ) 

holds where <p and <p' are, respectively, the perturbation po- 
tentials for the forward and reverse motions. The region of 
integration is determined by the area occupied in the xt 
plane by the wing, and the Kutta condition is assumed to 
apply to the trailing edge when in subsonic flight. 

If the left side of equation (19) is integrated by parts, the 
general relation becomes 

—^^w'dxdt=^*~-wdxdt (20) 

If the motion of a three-dimensional wing is to be studied, 

equations (19) and (20) must be modified to include an inte- 
gration with respect to y. 

Two further changes in equation ( 20 ) serve to simplify 
applications. In the first place, asymmetry is restored to 
the expression if two distinct systems of axes are used as in 
equation (14); in the second place, the pressure relation 

dip ^ /oi\ 

#o= — Po^7= — PoOo-gj (21) 


where p 0 denotes undisturbed pressure, permits the introduc- 
tion of pressure p in the integrands. The final expression 
for the three-dimensional case is, therefore, 


JJJptCzi’l/i, ti)w 2 (x h y u t l )dx 1 dy l dt 1 

= JJJ.fe, Vt, h)Wi(Pt, Vi, t^dx 1 dy 2 dt 2 (22) 


where the two motions now follow the same path in reverse 
directions but are referred to the two sets of oppositely 
oriented axes satisfying the relations 


*i— - £a+£, Vi = — 2/a+i?, 21 — Zj, t\ — — f2+F (23) 

where £, 77 , and r fix the relative positions of the two origins. 
Figure 3 indicates one possible orientation of the axes. 

Equation ( 22 ) reduces to a much simpler form, provided 
further restrictions are imposed on the upwash functions 


z, 



'Uh&u Vu U) and Wife, y it U). In order to fix the idea, con- 
sider the case in which the wings have traveled a time T', 
(T'—T/oq) and a distance U 0 T'=M 0 T. Let the two sys- 
tems of coordinate axes be placed such that t,=0 sets the 
starting time of the forward motion and t«=0 sets the start- 
ing time of the reverse motion; the two origins are further- 
more oriented such that they are at opposite ends of the 
root chord of the common plan form. Equation (22) then 
becomes 


dfij Jpi(xi,y u U)vh(x u y x , QdxAji 


PdU) 


=J o duj f 2Lfe, y», U)Wi fe, 2/2, ti)dxidyi 
' ° P s &) 


where the functions w t and vm have an implicit dependence 
upon T. If w 1 and vh remain constant for r 1 +ilfo^x= const, 
or const., the expression 


J JpiixuyuTtvisixuy^Tjdxidyi 

Pi{T) 

=J J Piixi, yi, T)Wi(x 3) yi, T)dx 3 dy 3 (24) 

PAT) 

follows after taking a derivative with respect to T of the 
original equality. Equations (24) and (14) are now equiva- 
lent in form, with T taking the role of an auxiliary param- 
eter. In this way, certain classes of unsteady motions can be 
treated simultaneously with steady motions. 

In the applications to follow it will be convenient to intro- 
duce into equation (24) upwash functions of the indicial 
type; that is, functions that are zero up to a fixed time and, 
after experiencing a finite discontinuity, remain constant 
for all subsequent values of time. Such indicial or , step 
variations can be assumed, say, for angle of attack, rate of 
pitch, and rate of roll since they satisfy the requirements 
underlying the derivation of equation (24). This choice of 
functions will prove to be advantageous in that the integrals 
of the responsive pressures will yield results relating the 
wing characteristics. Theorems to be given later will speak 
specifically of steady and indicial motions. It is to be under- 
stood, however, that the indicial results can be further 
extended when the same wing is assumed to be executing 


321000 — on 18 
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identical motions in forward and reverse flight. Thus, by- 
means of Duhamel’s integral (see, e. g., ref. 20), if f(t) is 
the response in the wing characteristic to a step variation 
in w at time t= 0, the response to an arbitrary variation 
with time of w can be written 


the product of the pressure and the slope in the x direction 
of the wing surface; when the wing has blunt edges with 
slopes having square-root singularities, these singularities 
yield an added contribution (ref. 19). In general, therefore, 
the drag D of a symmetrical section is given by 


F^J^t-r^dr (25) 

If it is known, for example, that' the lift per unit angle of 
attack is the same at corresponding values of time for a 
wing experiencing an indicial angle-of-attack change in 
forward and reverse flight, it follows that the build-up of 
lift is the same at corresponding values of time for all for- 
ward and reverse motions, provided the time histories of 
the motions are the same. The equivalence of lift would 
thus be established, for instance, for oscillatory variations 
in angle of attack. 

An alternative study of reverse-flow theorems for oscilla- 
tory motions could be based upon the modified wave equation 

Xcr+Xw+Xa+" 2 X=0 

which results from setting ^—e tat x in equation (18). Such a 
study would corroborate the conclusions drawn from equa- 
tions (24) and (25). 


APPLICATIONS 

The results of the foregoing analysis may be employed to 
determine a number of special theorems that are particularly 
useful in the calculation of the aerodynamic characteristics 
of twisted wings and of wings in nonuniform downwash fields. 
The theorems apply equally to wings acting either alone or, 
in certain cases, in combination with other wings or with 
cylindrical bodies having their generators alined with the x 
axis. Moreover, they apply not only to wings in steady 
motion but also to wings performing unsteady motions of the 
indicial type, or unsteady motions derivable therefrom. For 
wings in more complex unsteady motions, however, it will 
be necessary to refer to the more general results of equation 
(22). Some problems of this nature will be described at the 
end of this section. 

The applications to be included are exact within the frame- 
work of linear theory and involve no further rectrictions on 
the wing plan form or Mach number except in certain indi- 
cated cases where it will be convenient to use results based 
on slender-wing theory. The examples axe intended to be 
representative in nature. 

REVERSAL THEOREMS— STEADY AND INDICIAL MOTIONS 


D=D,+2 



dS 


(26) 


where D, is the drag attributable to the edges. 

If the subscripts 1 and 2 refer to the same wing in forward 
and reverse flow, respectively, and with the two systems of 
axes introduced in equation (23), local slopes are related as 
follows 

dz^XuyuT) 

, dxi 


dz?(xt,y 3 ,T) 


dx j 


(27) 


Equations (26) and (27), together with the reciprocal relation 
(24), yield 





- 2 JMS).^ 


(28) 

p* 

P, 


Since the geometry of the wing is fixed, the edge contribu- 

tions are the same, 


(29) 

and, consequently, 

Di—D 3 

(30) 


which confirms the relation stated in reference 9. 

THEOREM: The pressure drag in steady or indicial 
motion of symmetrical nonlifting wings is the same in 
forward and reverse flight. 

LIFT ON FLAT-PLATE WINGS 

The lift A of a wing may be determined by integrating the 
differential pressure Ap=p ( — ’over the wing plan form, 
thus 


For flat-plate wings, the local angle of attack of the wing sur- 
face is a constant 



ai (xi,yi,T)—ai= const. , a 2 (x 3 , 2 / 2 ,T)=a a = const. (32) 


Reversal theorems are defined here as relations between 
the aerodynamic characteristics of identical wings executing 
the same type of motions in forward and reverse flight. The 
results presented in this section apply not only to single 
wings in steady motion but also to combinations of wings, as 
in cascades or multiplanes, performing either steady motion 
or motions of the indicial type. 


Application of equations (31), and (32), and (24) yields the 
following: 


Tj\cc% — C C AjpiGC2 dSi — J^*Ap 2 ai(fiS , 2 — Is%<X\ 
Pi P 2 


or 


Iri/oCl — Aj/a 2 


(33) 


DRAG OF SYMMETRICAL NONLIFTING WINGS 

The drag of a symmetrical, sharp-edged wing in linear 
theory may be determined by integrating over the plan form 


THEOREM: The lift per unit angle of attack of flat- 
plate wings in steady or indicial motion is the same in 
forward and reverse flight. 
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This theorem generalizes the relation previously given by 
Brown (ref. 5) for steady motion. 


ai ~ it* 1 ’ aa ~U 0 Xa ~U 0 ^- 


(38) 


DAMPING IN ROLL OF FLAT-PLATE WINGS 

The rolling moment L' exerted on a wing, following the 
usual sign convention, is given by 


1 £'=— J* j'yApdS 


(34) 


The local angle of attack due to rotation about the x axis is 


’ 2&! 

1 U 0 ’ 


a-2- 


ViVi 




u 0 


Uo 


(35) 


where p’ is the angular velocity of roll, assumed constant. 
Application of equations (34), (35), and (24) yields the 
following: 
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THEOREM: The rolling moment per unit angular 
rolling velocity of flat-plate wings in steady or indicia! 
motion is the same in forward and reverse flight. 


DAMPING IN PITCH OF FLAT-PLATE WINGS 

Consider a wing, first, in forward flight and pitching with 
a uniform angular velocity (p about a lateral axis; second, 
in reverse flight and pitching with angular velocity q 3 about 
another lateral axis. Place each wing in a coordinate system 
such that the y axis coincides with the axis of rotation and 
designate the distances to the moment axes by Xo with 
proper subscripts, as shown in figure 4. In such a coordinate 
system, the pitching moment M Xq exerted on a wing, following 
the usual sign convention, is 


M Xq = — jj (x—x 0 )ApdS 


(37) 



Fioube L— Coordinate Bystems and symbols used In discussion of reverse-flow theorems 

for pitching wings. 


Application of equations (37), (38), and (24) yields the fol- 
lowing: 
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(39) 


This equation indicates that the pitching moment due to 
pitching velocity is, in general, not the same for wings in 
forward and reverse flight. However, if a^ 1 =x 0 2 =f, the 
pitching moment per unit angular pitching velocity of flat- 
plate wings in steady or indicial motion is invariant. 


SPECIAL RECIPROCAL THEOREMS AND APPLICATIONS 

In the following section, several special reciprocal theorems 
will be derived and applications will be illustrated. Recipro- 
cal theorems, in contrast to reversal theorems treated in the 
preceding section, are defined here as relations between the 
aerodynamic properties of wings in forward and reverse flight 
that have dissimilar camber, twist, and thickness distribu- 
tions but have the same plan forms. The motions may or 
may not be similar, although it is assumed in this section 
that both wings are in either steady motion or unsteady 
motion of the indicial type. As noted in the preceding sec- 
tion, the results apply equally to wings acting alone or in 
combination. 


SYMMETRIC NONLIFTING WINGS— STEADY MOTION 

The problems of paramount interest in the application of 
the general relations are found from considerations of pres- 
sure integrals over lifting surface; such problems will be 
given detailed treatment later. In the present section, a 
brief indication is given of the manner in which useful results 
can be derived for symmetric wings. The discussion will 
be limited to steady-state, two-dimensional, subsonic pres- 
sure fields although fairly obvious extensions can be carried 
out. 

If the geometry. of a real symmetric airfoil is prescribed 
the theoretical pressure distribution exists and is unique. 
If, however, the pressure distribution is prescribed, a real 
airfoil does not necessarily exist, but by means of reciprocal 
relations it is possible to derive certain conditions of com- 
patibility that need to be imposed. Consider, therefore, 
the two subsonic solutions 

Ui(x 1 )=0 > Wi{x\)=U ol-yfaF^x^) — c<a:i<o 


The local angle-of-attack distributions due to rotations 
about the y axis are 


and 


Vvixt), iipfe) 
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The first solution has square-root singularities in w at each 
end of the airfoil-and, correspondingly, singularities in u occur 
just ahead of the point x= —a and just behind x=a. On the 
other hand, equation (14) certainly applies if w* is zero at 
x=±a. If the origins of the two systems of axes are at 
the same position, it follows from equation (14) that ih(x) 
must satisfy the relation 

0= f (40) 

This result is useful in the calculation of airfoil shapes in- 
volving a change in pressure distribution from that of a known 
reference profile. The restriction on vh at the nose and tail 
implies that the derived and reference profiles have the same 
slope and radius of curvature at those points. The restric- 
tion on u- 2 , as given in equation (40), can be interpreted as a 
condition that must exist by virtue of the fact that the drag 
of an airfoil in two-dimensional potential flow in zero. 

As a second example, consider the solutions 

that represent velocity and slope of a thin ellipse of [thickness 
t and chord 2a. If toj is chosen as above, such that it van- 
ishes at the nose and tail of the airfoil, if is the correspond- 
ing velocity distribution, and if the two sets of axes are as 
before, equation (14) yields 

P“ dz * fa _ C a 
J dx I 1 J -a — xi 


Making the substitution Xi=— x% in this equation and inte- 
grating the left side by parts, one has 


£7oj^(2a;i— a.)z 2 J^ —2^ z^xjjdx^ 


P® 2x 2J rax 3 — ft 3 
j-a -yJaP—X 2 s 


u 3 (x 3 )dx 3 


For all real airfoils with cusped trailing edges, therefore, the 
area A 3 can be expressed as 



2x 3 -\-ax 3 — a 2 u 3 (x 3 ) 
ija 1 — x 3 U 0 


dx 3 


(43) 


LIFT — STEADY AND INDICIAL MOTION 

The reciprocal theorems offer considerable advantage in 
the calculation of the lift of wings having a nonuniform 
angle-of-attaek distribution or of wings in a stream having 
nonuniform flow directions. For these applications, it is 
convenient to consider a special form of the reciprocal theo- 
rem which relates the lift on a wing having arbitrary distri- 
bution of local angle of attack to that of the flat-plate wing 
of identical plan form in flight in the reverse direction. 
Since the solution of this latter problem is often known or 
• can be found relatively easily, the solution of the original 
problem is facilitated in' many instances. 

lift of arbitrarily cambered wings. — Consider two wings of 
identical plan form in flight in opposite directions, as shown 
in figure 5. Wing 1 is arbitrarily cambered and twisted and 
wing 2 is flat. 

<*i=cti(xi,yi,T), a 2 = const. (44) 


From this result, together with the general closure condition, 

/.‘.I (41 > 


a necessary condition for the closure of the second airfoil is 


P“ Vl 

J-a ^/a 2 — Xj* 


(42) 


As a final example, consider the solutions for —a<C.Xi<Ja 
, . 2tU 0 2 xi 2 —ax 1 —a i 

Wi fo) > 

3i/3a 2 -Ja 2 — x 2 


representing velocity and slope of a thin Joukowsky type 
airfoil. In this case, Wi vanishes at the tail and the down- 
wash distribution w? for the reverse wing may have a square- 
root singularity at the nose. The nose of the first wing is, 
however, blunt and for equation (14) to apply the second 
wing must have a cusped taiL Under these conditions, 
equation (14) yields 



2tUo_ r a 
3^/3 a 2 J-a 
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Figure 5.— Coordinate systems and symbols used In discussion of relation botwcon lift of 
arbitrarily cambered wings and loading on flat-plato wings. 


Application of equations (44) and (24) yields the following: 


or 


aj^i — \dS\ — i aiAp^dSi 


Pi 



(45) 


THEOREM: The lift in steady or indicial motion of a 
wing having arbitrary twist and camber is equal to the 
integral over the plan form of the product of the local angle 
of attack and the loading per unit angle of attack at the 
corresponding point of a flat-plate wing of identical plan 
form in flight in the reverse direction. 
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Equation (45) may be used to derive Munk’s integral for- 
mula for the lift of an arbitrarily curved airfoil in subsonic 
flow. Consider airfoils 1 and 2 placed in their respective 
coordinate systems, as indicated in figure 6. The angle-of- 
attack distributions on the two airfoils are given by 

dz i 

“ 1— dxi 

— a a — a a 

x \ x ” x z 

Fiqube G.-— Sketch of arbitrarily cambered airfoil Illustrating symbols used In equations (46) 

through (60). 

The loading per unit angle of attack on airfoil 2 is 


a 3 = const. 


(46) 



*2 



Figure 7.— Sketch of supersonlo-edged wing Illustrating symbols used In equations (51) 

and (52). 


sweep angle of the straightedge is A, the load distribution 
per unit angle of attack of the reversed wing is 


Apt _ 4g 0 cosA 
“2 cos 2 A 

and the lift is 


r. t = 4g -° — ff^ Idxtdy- 

Vsec 2 A — Mi JJ dx\ 


(51) 

(52) 


Ap2__4£o l a— r2 == 4go la+x 1 
a 2 ~ 1 3 \ a-\-x i p \ a—x i 


(47) 


whore g 0 I s free-stream dynamic pressure (^pqUq 2 ). Substi- 
tution into equation (47) yields the lift formula 


T _ 490 r° i a + x i 

Ll ~ p J-adx^a- 


Xi 


dxi 


(48) 


The corresponding formula for the lift of a tier of curved 
airfoils may also be derived similarly from the expression 
for the loading on an equivalent tier of flat airfoils. Con- 
sider, for example, an unstaggered lattice of flat-plate 
airfoils arranged vertically. If the gap distance between the 
plates is h, the loading per unit angle of attack is 


A p 2 4 q 0 ,/ira\ /sinh [ir(a—xj/ph] 
a 2 p VpfJ\ sinh Ma+zjfph] 

The formula for the lift on one of a lattice of identically 
cambered airfoils 2 is therefore 


4go ao . h /W\ f° dzi / sinh [vla+xdlPh] 

' l p \ph)j-a dx } \ sinh [*•(«— x{)/ph] ■ 


The load distribution per unit angle of attack for a two- 
dimensional supersonic wing is 


Ap 2== 4£o 
®2 P 


and, from equation (47), the lift is 



X 


a dzi 
-a dx 


(49) 


(50) 


The extension of this result to include supersonic-edged wings 
straight trailing edges leads to a result given originally by 
Lagerstrom and Van Dyke (ref. 21). If, as in figure 7, the 

* This result, as well as the detailed pressure distribution, has been derived by Mr. Paul 
F. Byrd of the Ames Aeronautical Laboratory by means of a direct Inversion of the singular 
Integral equation relating the aerodynamic pressure and geometry of the ftlrfoIL His work 
has also been used In deriving equations (80) and (81). 


A less obvious application yields the build-up of lift with 
time of an arbitrarily cambered two-dimensional supersonic 
airfoil starting impulsively from rest at a constant speed. 
Rewriting equation (45), lift is 



a-Mo T dz^T) Ap 2 (Xo,T) 
MoT dx i a 3 


dxi 


(53) 


Figure 8 presents sketches showing the final positions of the 
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Figdp.e §.— Regions In which different forms of equation ( 61 ) apply. 


airfoils relative to each other for various values of T. In 
reference 22 the expressions for the loading in equation (53) 
are given. Over the intervals denoted by ©, ®, and © in 
the sketch, these expressions are 


Region ® 
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lift on a wing in a nonuniform downwash field, — The 
reciprocal theorem of equation (45) can also provide a par- 
ticularly good method of determining the lift on a wing in 
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certain nonuniform downwash fields of known structure. 
Sucli problems arise whenever a wing acts in the presence 
of other wings, bodies, or propellers hut is always of prime 
concern in the determination of the lift on a tail acting in 
the downwash field of a wing. In most problems, the down- 
wash velocities at the position of the tail may be considered 
to be constant in the longitudinal direction and to vary in 
the spanwise direction, thus 


Substitution of equations (57) and (58) into equation (66) 
yields the following formula for the lift: 


T Li)- 

irU(jS a 2 
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(59) 


and 


ai=ai(yi,T), a 2 =const. 

Ll= ff ai 6?) S) dV ' 

Pa 


(55) 

(56) 


where la is the span load distribution associated with the 
load distribution A -p 2 . Summarizing, the lift in steady or 
indicial motion of a wing in a downwash field which varies 
across the span is equal to the integral over the span of the 
product of the local angle of attack and the span loading per 
unit angle of attack at the corresponding spanwise station 
of a flat-plate wing of identical plan form in flight in the 
reverse direction. This statement generalizes the result 
given recently by Alden and Schindel (ref. 23) for steady 
flow about wings having supersonic leading and trailing 
edges and stream wise ^ido edges. 

As for example, consider the problem of determining the 
lift on a wing at a geometrical angle of attack of zero resulting 
from the presence of an infinite line vortex of strength r 
extending in the flight direction. The wing will be considered 
to have such a plan form that its span loading when in 
flight in the reverse direction is elliptic. The notation is as 
shown in figure 9. For this problem, therefore, the span 



loading of the wing in reverse flight is given by 


la 2Lj I ^ « Va 

a-a tt 8 a 2 V S 2 


The local angle of attack of the original wing due to the 
presence of the vortex is given by 


<*i= 


r - yi-p 
’2-rI7 0 l(yi-py+'r 


The lift on a wing in the vicinity of a number of such 
vortices may be found by superposition. 

lift due to deflection of a portion of the wing surfaoe. — 
Let a portion P' of the surface of wing 1 be deflected a con- 
stant angle 5 and the remainder of the wing, be a flat plate 
alined with the free-stream direction. Let wing 2 be a 
flat-plate wing inclined at an angle of attack a, thus 


(5 on P' 

cci=< a%— const. 

(0 elsewhere 


(60) 


Substitution of equation (60) into (45) yields 
result: 



the following 

(61) 


The lift in steady or indicial motion per unit angular 
deflection of a portion of the wing surface is thus equal to 
the lift per unit angle of attack on the corresponding portion 
of a flat-plate wing in flight in the reverse direction. This 
generalizes a result given previously by Morikawa and 
Puckett (ref. 24) for steady flow about low-aspect-ratio wings. 

This rule is very useful in the determination of the lift 
resulting from the deflection of a flap or control surface. 
This is particularly true for supersonic speeds since the load- 
ing on the related flat wing is often a constant over a largo 
portion of the area of integration. 

As a further example, consider the case of a low-aspect- 
ratio wing having a straight trailing edge and mounted on 
an infinite cylindrical body of revolution. The entire wing- 
body combination is at zero angle of attack except for the 
flaps on the rear of the wing that are deflected an angle 5. The 
problem is to determine the lift on the entire wing-body 
combination due to the deflection of the flaps. Slondor- 
wing-theory results of reference 25 are to bo used. The 
notation is indicated in figure 10. 



(58) 


Figure 10. — View of lifting slender wing-body combinations. 
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The solution of this problem is particularly facilitated by 
the fact that slender-wing theory indicates that the loading 
on wing 2 is concentrated on the leading edge, as shown in 
the figure. Therefore, the lift of wing 1 is found by inte- 
grating the span loading curve of wing 2 over the portion of 
the span between y=j and y=s. Thus 

<62) 

flap 

The span loading on wing 2 is given by 

k (Hi) =42o«iO-^(t-Q (S~$) (63) 

The lift due to the deflection of the flap is therefore 


X,,.= 
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a 1 

-5 arc sin 
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1 — (a/s)* 

[l+(am<J/s) 2 - 2 (a/sy 

(fls) 2 [l-(am 


! 


( 64 ) 


A plot of the results is shown in figure 11. The lift per unit 
angle of flap deflection (A/ 5 ) /Iap has been nondimensional- 
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Figure 11.— Lift of a slender wing-body combination resulting from flap deflection. 


ized by dividing by the lift per unit angle of attack (L/a) w - B 
of a slender wing-body combination of identical plan form. 
From reference 26 , (L/a)jy- B is given by 

(Zr/a) I jr_ B =27rg , o^ a (^ 1 — ^ 5 ^ (65) 

ROLLING MOMENT— STEADY AND INDICIAL MOTION 

The calculation of the rolling-moment characteristics of 
wings having a nonuniform angle-of-attack distribution or 


of wings in a stream with nonuniform-flow directions can be 
performed in many cases through use of the reciprocal 
theorem in a manner analogous to that described for the 
lift characteristics in the preceding section. In every case, 
the rolling moment of the given wing will be related to the 
lift on a rolling flat-plate wing. For the sake of simplicity, 
all the present examples will be confined to the case where 
the rolling moments are evaluated about the x axis, con- 
sidered to he in the plane of symmetry. 

Rolling moment of arbitrarily cambered wings. — Consider 
two wings of identical plan form in flight in opposite direc- 
tions, as shown in figure 12. If the local angle-of-attack 



distribution of wing 1 is arbitrary and that of wing 2 varies 
linearly with y (which might be likened to either a wing with 
linear twist distribution or to a flat-plate wing rolling about 
the x axis with constant angular velocity p/), 

aci=ai(x u yi,T), a 2 =fl^~ C66) 


and the following relations can be written: 

Ifl = — JJ* yi&Pid3i.=jJy-ApflSi == JJ* P? a ^PflSi= 

Pi P 3 


(67) 


THEOREM: The rolling moment in steady or indicial 
motion of a wing having arbitrary twist and camber is 
equal to the integral over the plan form of the product of 
the local angle of attack and the loading per unit (pi'/U 0 ) 
at the corresponding point of a rolling flat-plate ■wing of 
identical plan form in flight in the reverse direction. 

Applications of this theorem follow in a manner very 
similar to that described previously for the- corresponding 
theorem regarding lift. 

Rolling moment on a wing in a nonuniform downwash 
field. — Consider a wing placed in a flow field in which the 
downwash velocities at the position of the wing are constant 
in the longitudinal direction and vary in the spanwise direc- 
tion. The related wing is again a flat-plate wing rolling with 
an angular velocity p 2 ' as described in the preceding section, 
thus 

<*i=“i(2/i,T), (68) 
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The rolling moment of the first wing is then given by 

i/== f /“ 2 (pTm) ds ^ /-, (jJjud dy * (69) 

P, 


or, in words, the rolling moment in steady or indicial motion 
of a wing in a downwash field which varies across the span is 
equal to the integral over the span of the product of the local 
angle of attack and the span loading per unit (p//Z7 0 ) at the 
corresponding spanwise station of a rolling flat-plate wing of 
identical plan form in flight in the reverse direction. 

Eolling moment due to deflection of a portion of the wing 
surface. — Let a portion P' of the surface of a wing be de- 
flected a constant angle 8 and the remainder of the wing be a 
flat-plate alined with the free-stream direction. The related 
wing is a flat-plate wing rolling with angular velocity p/ 


«i— 


5 on P' 

0 elsewhere 




pi Vi 

" u 0 


(70) 


Substitution from equation (70) into (67) yields' the following 
result: 



(71) 


Thus, the rolling moment in steady or indicial motion due 
to a given angular deflection of a portion of the wing surface is 
equal to the lift per unit (p/l U 0 ) on the corresponding portion 
of a rolling flat-plate wing of identical plan form in flight in 
the reverse direction. 

As an example, consider a wing-body combination con- 
sisting of a low-aspect-ratio wing having a straight trailing 
edge mounted on an infinite circular cylinder, as shown in 
figure 13. The body is at zero angle of attack, the right wing 


•?2 



Figube 13, — Sketch Illustrating symbols used In discussion of rolling moment resulting 
from differential deflection of wings of slender wing-body combination. 


P' is deflected an angle 8, and the left wing P" is deflected — 5. 
The problem is to determine by means of slender-wing theory 
the rolling moment exerted on the entire wing-body com- 
bination. The notation is indicated in figure 13. 


{ —5 on left wing, P" 
0 on body 

+5 on right wing, P' 


a 3 


U 0 


(72) 


wing 1 can be found by integrating the span loading on wing 

2. 


(^fm) dS >— s f_, (^ju) d y>+ 

p% 

s Slim dy '= 2s rXim) dv ‘ 


where, from reference 26, 
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fr'IUo 


«sx{(l+; arc cos 
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vV yJ (^-^(s’+^J 

(73) 

The resulting expression for the rolling moment is 


V=4«„S^{l^,i)[^(l-^+lr)(l+§orocosj^)^ 

lr “ -*>]+*<* *> [I (h^+fX 1 + 

! arc cos j^)+£ S(1 -S’jJ+f? jq^+ 
^a-JP)(l+|arcco, T ^)J (74) 

I 

where 

R=a/s 

. ■ 1 

p=arc sin - 7= ===== 

XT+IP 

A plot of the results is shown in figure 14. The rolling 
moment has been nondimensionalized by dividing by the 
value corresponding to that of the wing alone (R— 0). 


PITCHING MOMENT— STEAD V AND INDICIAL MOTION 

A number of useful relations regarding the pitching- 
moment characteristics of wings may be found through 
application of the reciprocal theorem. Since the general 
procedure is closely analogous to that of the preceding sec- 
tions, the following 'discussion will be brief. 

Pitching moment of arbitrarily cambered wing. — Consider 
the problem of determining the pitching moment Mi about 
the origin of wing 1 possessing an arbitrary distribution of 
camber. The related wing in flight in the reverse direction, 
wing 2, is a flat-plate wing of identical plan form pitching 
about the moment axis of wing 1, as indicated in figure 16, 
thus 

q 3 Xi _q 3 ( x 3 —£) 

Uo U 0 


Since slender-wing theory indicates that the loading on wing 
2 is concentrated on the leading edge, the ro lling moment of 


ai=ai(xi,yuT), 


a 2 = 


(75) 
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Fioobe 14.— Rolling moment of slender wing-body combination resulting from differential 

deflection of wings. 


The necessity for pitching wing 2 about the moment axis 
of wing 1 may be removed by considering wing 2 to be 
re-expressed in terms of two component wings having angle- 
of-attack distributions given by 

qAxr— Xo»>) qAx<n — I) * 

oy=- 77 > a 2 "= 77 = const. (77) 

Uo U 0 


Wing 2 ' is thus pitching with angular velocity g 2 > about an 
axis at rs'= 2 ta' and wing 2" is a flat-plate wing at a con- 
stant angle of attack. The pitching moment on wing 1 is 
then given by 


Pi' Pi" 


“If- (| w) -»*■+(*--« >If« (^r) ds >" (7S) 

Pi- Pi" 

Applications of pitching-moment theorem. — The applica- 
tion of equation (76) or (78) to problems analogous to those 
discussed in the preceding section can be carried out in a 
straightforward manner. Consider, first, unstaggered lat- 
tices on airfoils such that the airfoils in lattice 1 have arbi- 
trary camber distributions and those in lattice 2 are flat 
plates pitching about their midchord positions. The angles 
of attack in the two lattice systems are 


a 2 (Xi,Zi)~22*l ( 79) 



The pitching moment of wing 1 is given by 

u ‘— !J^' dS '=J!Awir) dS ' 

SS a '(-m) iS - m 

P 2 


and the load distribution on each airfoil in lattice 2 is, in 
subsonic steady flow, 

Api(xJ=q 0 (' arc cos J^sech cosh ( 80 ) 

where 2a is chord length. Equation (76) yields, for pitch- 
ing moment of the first airfoil about its midchord point, 
the result 


M ‘ cos [ sech (M) “•‘‘fir)]*’' (81) 


A second example, illustrating unsteady effects, is the fol- 
lowing: Let wing 1 be a flat-plate wing, then «i is constant, 
and equation (78) simplifies to 


Lfi 

1 qr/Uo 


OC2" 


(82) 


where Li is the lift on wing 2' pitching about xP =x-A , and 
Li" is the lift on an inclined flat-plate wing. Equation (82) 
may be expressed in terms of conventional stability deriva- 
tives as follows: 


THEOREM: The pitching moment in steady or indicia! 
motion of a wing having arbitrary twist and camber is 
equal to the integral over the plan form of the product 
of the local angle of attack and the loading per unit 
(qz/Uo) at the corresponding point of a flat-plate wing 
of identical plan form in flight in the reverse direction 
and pitching about the moment axis of the first wing. 


(C ma \=(C L! )i'+(^-^yC La )i- (83) 

An application of this result to unsteady-flow problems is 
indicated in figure 16 obtained from indicial-lift and pitching- 
moment results of reference 27. This figure shows the 
growth of lift and pitching moment on triangular wings with 




REPORT 1119 — NATIONAL ADVISORY COMMITTEE EOR AERONAUTICS 


266 


2.5 


2.0 


1.5 


-.5 


- 1 . 0 . 




M 0 m Z 

' 


M C L a )z m 

■ 







. 

— "" 

*■ 

' i C l-q)z' 




&o)| : 










’ 




• 











.6 


2T0 


25 


1.0 1.5 

Chord lengths traveled 

FiOTOE IB. — Indtclal lilt and pitching moment on triangular wings with supersonic edges. 


supersonic edges at a Mach number of 2 following indicia! 
angle-of-attack and pitching-velocity changes. In these 
results, the rotation and moment axes are always at the 
leading edge or apex, therefore, x^'—Q and £=Co. It may 
be seen that the three curves are related in the simple linear 
manner indicated by equation (83). 

If ai is independent of x and varies only in the spanwise 
direction, that is, if ai=ai(yO> the pitching moment on 
wing 1 is given by the following equation, analogous to 
equation (56) for lift: 


Mi= /_, ai (}jtm) (iO (84) 


If a portion P' of the surface of wing 1 is deflected a con- 
stant angle S and the remainder of the wing is a flat plate 
alined with the free-stream direction, the following relations 
hold: 

|5onP' 

«i=< 

( 0 elsewhere 


and 


P'r Pv ' 


( 85 ) 


RECIPROCAL RELATIONS INVOLVING MOTION INTO A GUST 

All previous applications that have been considered wore 
derived from equation (24). In the present section, the 
more general equation (22) will he used to develop two the- 
orems which relate the build-up of lift on a wing entering 
a gust and the build-up of circulation on the same wing 
moving indicially but in the opposite direction. The rela- 
tions to be obtained hold for the Mach number range for 
which the wave equation applies. Under the special assump- 
tions of incompressible flow, the results in two dimensions 
establish a direct connection between the circulation function 
calculated by Wagner (ref. 28) and the gust lift curve cal- 
culated by Kussner (ref. 29) and von Kfirmfin and Sears 
(ref. 30). A proof of the connection between these functions 
for two-dimensional incompressible flow has been given by 
Sears (ref. 31). 


TWO-DIMENSIONAL FLOW 

A flat plate is assumed to be moving in two modes of 
motion: In the motion associated with the x^,Zi,tq, axes, the 
wing starts at' time zero (t»=0) and moves at a constant 
velocity U 0 and at a constant angle of attack; the motion 
associated with the Xi,z u t u axes starts at time zero (h=0) 
with the wing moving'in the opposite direction at a velocity 
Z7 0 and entering a sharp-edged gust. The gust exists for 
all Ti less than zero and has a vertical velocity w t = — a s U 0 . 
The two wings, therefore, have angles of attack as follows: 

a 2 {xt, ^)=a 3 =const. for— Afo<a<a^<2a— M 0 h, 0 

f— A/ o h<ri<0 0 ^ h g 2a/M 0 

aifo, <i)=«*=const. for] 

I— M 0 ti<xi<2a— M 0 ti 2a/M 0 ^ti 

The two-dimensional form of equation (22) yields 

“ 3 Jo dk L-Z^ A ^ dS 

A (86) 

where the region A is bounded by the lines 21=0, Ti=2a— 
Mdi, ti=T, and x x = The integral on the right can be 

rewritten as a line integral by means of the identity 

— J*P cos ( t , n)ds=jj^~ dS 

and equation (86) becomes 


a 3 f Lg (ti)dti = po^o^r J ~ Atpj cos (fj, n)ds 
O 


where the line integral extends around the boundary of the 
region A. Since A^a vanishes on the lines x i =—M 0 ti and 
f 2 = 0 , the equation becomes 


a=a 


Jo I'ziWti-P'Ao** J 2a _ MoT ^( x *’ 


Differentiation with respect to T yields 


a 1 L t (T)=p 0 Uooc t A<p 2 (2a—MoT, T) 
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The discontinuity in <p 3 is evaluated at the trailing edge at 
time T and is therefore equal to the circulation r 2 of the air- 
foil at time T. The equality thus becomes 


Ml 


PoUo 


UT) 

a 3 


(87) 


THEOREM: The circulation per unit angle of attack 
of a flat plate moving indicially with a velocity U 0 is 
proportional to the lift per unit a t of the plate entering 
a sharp-edged gust having a uniform vertical velocity 
equal to w s =—a 1 U 0 . 

In figure 17, the time variation of these variables, as well 
as the lift of the indicia! wing, is indicated for low speed and 
for flight Mach numbers equal to 0.8 and 1.46 as determined 
from references 22 and 27. 



Half-chords traveled Half-chords traveled Half-chords traveled 
M 0 M 0 M 0 *O.Q ' M 0 * 1.46 

Figure 17. — Growth of ci a and cu with chord longthg tiavoled. 


THREE-DIMENSIONAL FLOW 

The extension of the above results to three dimensions 
follows directly. The origin of the x u y h z x , t v axes is assumed 
to be initially at the foremost point of the wing in reverse 
motion. The two wings have, respectively, angles of attack 
a 2 =const. over the reverse moving plan form for all values 
of, time and ag=const. over the region occupied simultane- 
ously by the forward moving wing and the gust. Equation 
(22) gives 


wing be given by the equation 

Xi=f(±yi)—M 0 ti or yi= ±s(x 1 +M a t 1 ) . 

For the reverse wing and its coordinate system, this edge, 
which is now the trailing edge, is 

x 2 =c T —M<h—f(±y i ) or y J =±s(c T — a*— M<&) 

where c r is the root chord. The reverse-flow integrals of 
equation (88) then become 

rr ce, rs(Cr~xj 

a i L/ti)dti= podoctg drj 

J0 JCr-M 0 T 

A<P* [x s , I f 3 , Cr ~ X2 tif ±Vl) ] dy 3 

where s(x j) is the local half-span of the wing. Differentia- 
tion with respect to T yields 

aaig(T)=Pof7 0 a r J' A^? \f T ~ M 0 T, y%,T— dy 3 

(89) 

The discontinuity in <p is thus to he integrated spanwise at 
the rearmost point of the indicial wing; this follows from the 
relation 

A Pa [x,y,T-^g2)]=A^(x,2/,t); T-&^<t<T 

which fixes the vorticity in the wake of the wing once it is 
shed from the trailing edge. 

It remains to mention the nature of the limits ±s(M 0 T). 
As shown in figure 18, the. span width of the vortex wake at 
the trailing edge is, during the early stages of the motion, 
dependent on the local span width of the wing. The width 
2s(M 0 T) of wake is, in fact, equal to the maximum width of 
the portion of the first wing that lies within the gust. After 


jo ^ ll jj Ap ^ Xl ’ Vl ’ t ^dxidy 1 — 

m 


Pododg 


JIf- 


dpafoa, Vi, t») 
btj 


dxidy 3 dti 


( 88 ) 


The integral on the right can be rewritten as a two-dimen- 
sional surface integral by means of the identity 

— jj P cos (t, n)dS ^ dxdydt 

and equation (88) becomes 


0 : 2 ^ Lg(ti)dti — pQ<Xoa t A ^?2 cos (V 2 , Ti)dS 2 


where the integral on the right extends over the boundary of 
the volume in x?, y % , t 2 space occupied by the wing and the 
gust. The value of A <p 3 must, of course, vanish on the lead- 
ing edge of the wing and at < 2 = 0 . In order to fix the limits 
of integration, suppose the wing is symmetrical about its 
longitudinal axis and let the leading edge of the forward 



Figure 18. — Sketch illustrating nature of Integration limits In equation (89). 
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sufficient time has passed for the vortex wake of the indieial 
wing to- develop its full span width at the trailing edge, 
s(M 0 T) becomes s 0 or semispan of the wing. From equation 
(89) one may conclude the following: 

THEOREM : The lift per unit a v of a flat-plate wing en- 
tering a sharp-edged gust having a uniform vertical 
velocity equal to w t = — is proportional, at each 
instant of time, to the spanwise integral at the trailing 
edge of the vorticity shed -by the same wing moving 
indicially in the reverse direction with a velocity U 0 . 

As a direct example, this theorem has been used to confirm, 
from a knowledge of the indieial solution, the sharp-edged- 
gust lift of the rectangular-plan-form supersonic wing given 
by Miles in reference 32. 

Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Feb. 19, 1952. 
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